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Description 

Technical Field 

5 The invention relates generally to preserving the frequency information of a sampled analog signal and, more 

particularly, to a method and apparatus for interpolating between data samples such that the interpolated values contain 
the same frequency domain information contained in the sampled analog signal. 

Background of the invention 

10 

Analog signals contain frequency domain information and time domain Information that may be of interest to per- 
sons analyzing the signals. For various reasons, analog signals are often sampled by a time-based sampling process 
prior to analyzing the signal. Data samples produced by the sampling process contain both time domain information 
and frequency domain information about the analog signal. The time domain information, such as amplitude, for ex- 
15 ample, can be analyzed by measuring the amplitude of the data samples with an appropriate measuring device, such 
as a voltmeter, for example. Likewise, the frequency domain information can be analyzed by an appropriate apparatus, 
such as a digital spectrum analyzer, for example. 

In order to find time and frequency domain information for points between data samples, an interpolation technique 
must be used. The prior art offers numerous time-based methods of interpolation that provide varying degrees of 
20 accuracy between the interpolated values and the analog signal. Like the original data samples, time domain informa- 
tion contained in an interpolated value, such as amplitude, can then be measured by an appropriate device, such as 
a voltmeter. 

Likewise, if frequency domain information is sought for points between data samples, the data samples must be 
interpolated to find these values. Unfortunately, a time-based sampling process that produces data samples of an 

2S analog signal also produces spectral images of the analog signal's frequency spectrum. As a result, the data samples 
contain frequency domain information that includes the analog signal frequency spectrum information plus image in- 
formation. The image information corrupts the frequency spectrum Information so that the frequency domain information 
contained in the data samples does not accurately replicate the frequency spectrum of the analog signal. Prior art 
interpolation techniques further corrupt the frequency information contained in the data samples. As a result, the fre- 

30 quency information contained in the interpolated values does not accurately represent the frequency Information con- 
tained in the analog signal. 

In Schafer et al. "A Digital Signal Processing Approach to Interpolation," Proceedings of the IEEE. Vol. 61 , No. 6, 
pages 692-702 (June 1973), the authors describe a prior art interpolation system for resampling an input signal at a 
uniformly spaced resample times. The process involves convolving the input samples with a filter impulse response 
3S at the resample times and summing. Since the input and output sample rates are related by a constant, the value of 
the filter impulse response need only be known at a predetermined number of uniformly spaced points. This system 
does not allow interpolating at arbitrary resampling times, and it ignores the effect of the interpolation in the frequency 
domain. 

EP 099. 600 discloses a related system which provides an output signal sampled at an irrational multiple of the 
<o input signal sample rate. Again, this reference does not allow interpolating at arbitrary respamling times and ignores 
the effect of the interpolating operation on the frequency domain. 

As can be appreciated from the foregoing discussion, there is a need to provide a method and apparatus for 
interpolating between data samples such that the interpolated values contain the same frequency information contained 
in the sampled analog signal, i.e., such that the frequency spectrum of the sampled analog signal is presen/ed. The 
<s present invention is directed to a method and apparatus for interpolating between data samples using a digital filter in 
a novel manner that achieves this result, i.e., produces interpolated values that preserve the frequency spectrum of 
an analog signal. 

Summan/ of the Invention 

so 

In accordance with this invention, a method and apparatus for producing interpolated values that contain accurate 
information about the original frequency spectrum of a sampled analog signal is provided, as defined in claims 1 and 
5. respectively. 

As will be appreciated from the foregoing summary, the present invention provides a method and apparatus for 
55 presenting the frequency spectnjm of a sampled analog signal by interpolating data samples using a digital filter such 
that the interpolated values contain frequency domain information that can be used to accurately replicate the frequency 
spectrum of the analog signal. 
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Brief Description of the Drawings 

The foregoing, and other features and advantages of this invention, will become nnore readily appreciated as the 
same becomes further understood by reference to the following detailed description when taken in conjunction with 
s the accompanying drawings, wherein: 

FIGURE 1 is a simplified block diagram of a preferred embodiment of an interpolating apparatus formed in accord- 
ance with the Invention; 

FIGURE 2 is a simplified flowchart illustrating the functional steps of a program for controlling the apparatus ilius- 
10 tratedin FIGURE 1; 

FIGURE 3 is a simplified block diagram depicting an alternative embodiment of the apparatus illustrated in FIGURE 
1. 

FIGURES 4a, 4b and 4c are more detailed flowcharts of the functional steps illustrated in FIGURE 2; 
FIGURE 5 illustrates the relationship between frequency domain information contained in data samples and the 
15 frequency response of a suitably designed digital finite impulse response ("FIR") filter suitable for use in the ap- 

paratus illustrated in FIGURE 1 ; 

FIGURE 6, lines A-E, illustrates the interrelationship between the frequency response of a digital FIR filter and the 
continuous impulse response of the digital FIR filter; and, 

FIGURE 7, lines A-E, is a series of waveforms illustrating the convolution of the data samples with the continuous 
20 impulse response of a four-point digital FIR filter suitable for use in the apparatus illustrated in FIGURE 1. 

Description of the Preferred Embodiment 

There are situations where it is desirable to analyze, in the frequency domain, data that has been created in the 

25 time domain. Such a situation arises when an original analog signal is sampled by a time-based sampling process that 
produces data samples. While there are numerous interpolation techniques that will provide interpolated values in the 
time domain based on the data samples, in general, such techniques ignore the frequency domain information con- 
tained in the data samples. Asa result, the interpolated values may not contain accurate frequency domain information. 
As a consequence, there has developed a need for an interpolating technique that will produce interpolated values 

30 containing accurate frequency domain information. Among other benefits, such values can be used to analyze the 
frequency spectrum of the original analog signal. Another example of frequency domain information that is of interest 
to some industries is the changes in phase angle of the shaft of a rotating machine that occur as the machine is 
operated. The frequency spectrum of analog signals created when such a machine is operated contains shaft phase 
angle information. As will become obvious from the following discussion, the present invention provides a method and 

35 apparatus for preserving the frequency spectrum of an analog signal by interpolating data samples using a digital filter 
in a novel manner so as to overcome the foregoing problems of the prior art. 

FIGURE 1 Illustrates a preferred embodiment of an apparatus for interpolating data samples in accordance with 
the present invention. The apparatus illustrated in FIGURE 1 comprises: a signal processor 10; a sampler 12; a storage 
device 14; and a clock 16. The sampler 12, the storage device 14 and the clock 16, per se. do not form part of the 

40 present invention. They are included so that the invention will be more easily understood. An original fluctuating analog 
signal, designated Vq, which has an original frequency spectrum, is applied to the sampler 1 2. The sampler 1 2 samples 
the Vq signal and produces a digital data signal comprising data samples, designated Vg, taken at known sample 
times, designated Xq, The ts times are provided to the sampler 12 by the clock 16. The Vs values and the ts times are 
stored in the storage device 1 4. 

45 As will be better understood from the following discussion, the signal processor 10 receives control infonnation 

from an external source, such as time Inten/als between the ts times, designated At; resample times at which interpo- 
lation is to be performed, designated tf,; and. the number of tp times required, i.e., a data block size, designated Bis. 
As required, the signal processor 10 also receives the Vs and ts values stored in the storage device 14. Based upon 
the externally generated control information and the Vs and ts values, the signal processor 10 produces a series of 

so interpolated values, designated Viy. The V|v values contain frequency domain information as well as time domain 
information. As will be better understood from the following discussion, the V|v values can be used to accurately rep- 
licate both the time and frequency domain information contained in the Vq signal. The signal processor 10 can take 
the form of a specialized signal processor, i.e., a signal processor containing a program dedicated to carrying out the 
Invention, or a suitably programmed general purpose computer. 

55 The broad functional steps of the program that controls the signal processor 10 are illustrated in FIGURE 2 and 

discussed next. Initially, the signal processor 10 reads into its internal memory the Vs values, ts times, tpj times, At 
intervals, and the data block size, Bis. Next, the program causes the signal processor 10 to read the first t^ time and 
find a set of Vg values occurring at tg times that are located near the first tR time. As will be better understood from the 
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discussion below, the size of the set of Vs values is determined by the nature of a digital filter, specifically a digital finite 
impulse response (FIR) filter that, preferably, is implemented in software form in the signal processor 10. The digital 
FIR filter has a time domain or continuous impulse response function, designated h(x), where x is a normalized time 
(i.e.. X = t/T, where T is defined below and t is some time of interest, such as t^, for example) and a frequency domain 

s or frequency response function, designated H(y). where y is a normalized frequency (i.e., y = fT). The continuous 
impulse response function. h(x), has a finite width, T The size of the set of Vg values is determined by the impulse 
response width, T For example, if the impulse response width, T, encompasses four ts times, the signal processor 10 
will find the four Vq values associated with the four ts times nearest the first resample time. tf,. In this case, the digital 
FIR filter is commonly referred to as a four-point filter. 

10 After the set of Vs values has been found, the program instructs the signal processor to convolve the continuous 
impulse response function. h(x), of the digital FIR filter with each of the Vg values in the set. to produce a convolved 
value, designated V^. located at tp. for each of the Vg values in the set. As will become better understood from the 
following discussion, during this convolution step, the frequency response function, H(y). of the digital FIR filter performs 
a filtering operation. That is, the filter reduces the effects of spectral images introduced by the sampling process. Next, 

'5 the signal processor 10 sums the Vq values to produce a first interpolated value, V|v, i.e., the V|v value associated 
with the first tp time. Once this V,v value has been computed, the program instructs the signal processor 10 to read 
the next tp time, and the functional steps discussed above are repeated. The functional steps illustrated in FIGURE 2 
are repeated until the signal processor 10 has produced V,v values for each of the tp times received and read by the 
signal processor 10. As a result of the filtering operation discussed above and other aspects of the filtering operation 

20 which are discussed below, the V,y values contain accurate frequency spectrum information of the original signal (i.e., 

While the presently preferred embodiment of the invention comprises a signal processor 10 and a program for 
carrying out the invention, as discussed above and illustrated in FIGURES 1 and 2, it is to be understood that alternative 
embodiments exist which are also suitable for carrying out the invention. One such alternative embodiment is a hard- 

2S ware embodiment of the sort illustrated in FIGURE 3. In this alternative hardware embodiment, the sampler 12. the 
storage device 14 and the clock 16 perform the functions discussed above. The alternative apparatus comprises: a 
receiver 18; a data locator 20; a digital filter 21; a convolver 22; and, a summer 23. The digital filter 21 is preferably a 
digital FIR filter having the h(x) and H(y) functions noted above. 

The receiver 18 receives the Vg values and tg times from the storage device 1 4. The data block size. Bis, the time 

30 interval between tg times. At, and the resample times, tp, are provided to the receiver 18 by an external source. The 
data locator 20 finds the set of tg times nearest the first tp time and determines the associated Vg values. As noted 
above, the size of the Vg value set is determined by the impulse response width, T, of the impulse response function, 
h(x). of the digital FIR filter 21 . Once the Vg values have been found, the convolver 22 convolves the continuous impulse 
response function, h(x), with the Vg values to produce a Vc value located at the tp time for each tg time in the set. 

35 Next, the summer 23 sums the Vq values at the tp time to produce the V,v value associated with the first tp time. The 
apparatus illustrated in FIGURE 3 repeats the above process until V,y values have been produced for each tp time 
received by the receiver 18. In a manner similar to that discussed above for the signal processor 10 (FIGURE 1), the 
frequency response function, H(y), of the digital FIR filter 21 performs a filtering operation. As a result, the original 
frequency spectrum of the Vq signal is accurately reproduced in the frequency domain information contained in the 

40 V|v values. 

Before returning to the preferred embodiment of the invention (FIGURE 1) and discussing, in more detail, the 
functional steps of the program illustrated in FIGURE 2, a general 'background' discussion of various aspects of the 
present invention is presented so that the invention may be more easily understood. More specifically, the novel con- 
cepts of: using the continuous impulse response function, h(x), of a digital FIR filter to perform time-based interpolation 

45 of the data samples, Vg; and, using the frequency response function, H(y), of the digital FIR filter to produce V|v values 
that contain accurate frequency spectrum information of the original signal (i.e.. Vq). 

As is well known in the signal sampling art, digitally sampling an analog signal produces spectral images of the 
analog signal's frequency spectrum. As a result, the data samples produced by a digital sampling process contain 
frequency domain information relative to both the frequency spectrum of the analog signal and images of the frequency 

so spectrum. When the data samples containing such information are interpolated, a portion of the image information is 
contributed to the analog signal frequency spectrum information contained in the interpolated values. This contribution 
of image information results in aliasing errors that distort the frequency spectrum information contained in the interpo- 
lated values. As a result, the frequency spectrum information contained in the interpolated values cannot be used to 
accurately determine the frequency spectrum of the analog signal. As will become better understood from the following 

55 discussion, the digital FIR filter used in the present invention reduces the contribution of the spectral image information 
in the interpolated values. As a result of this and other filtering operations performed by the filter, which are discussed 
below, the frequency domain information contained in the interpolated values can be used to accurately detenmine the 
frequency spectrum of the analog signal. 
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FIGURE 5 illustrates the relationship between a digital FIR filter suitable for use with the present invention and the 
frequency spectrum of an analog signal and Its spectral Images. An original fluctuating analog signal, such as Vq, for 
example, may be described as having a frequency spectrum 24 with the general shape depicted in FIGURE 5. So as 
to simplify the following discussion, the frequency spectrum 24 is centered about a zero frequency point. As noted 
5 above, spectral Images are produced during sampling of the Vq signal. As was also noted above, frequency domain 
information relative to the frequency spectrum 24 and the spectral images is contained in the resulting data samples, 
Vg. As is well known in the sampling art, by sampling the Vq signal at a sampling rate, designated fg, using an over- 
sampling factor, designated m, where m > 1, the spectral images are centered at integer multiples of f^ {i.e., ilfg, ±2fg. 
±3fs. etc.). A representative pair of spectral images 25a and 25b are depicted in FIGURE 5 as dashed cun^es and are 
10 centered at ±1f^. As is also well known in the sampling art, the bandwidths of the frequency spectrum 24 and the 
spectral images 25a and 25b are the same and are a function of the sampling rate, fg, and the oversampling factor, m. 
Specifically, the bandwidths are equal to the sampling rate divided by the oversampling factor (i.e., f^/m). Thus, as 
illustrated in FIGURE 5, the frequency spectrum 24 has a bandwidth that lies between ±iJ2m. Likewise, the bandwidth 
of the representative spectral image 25a lies between + fg (^^) and + fg (^g^), and the bandwidth of the other rep- 



15 resentative spectral image 25b lies between -fg (^^) and -fg (%^)- As noted above, and as will be better understood 
from the following discussion, a filter, such as a suitably designed digital FIR filter, can reduce the contributions of the 
spectral images 25a and 25b to the frequency spectrum 24 when the interpolation of the V3 data samples is performed. 

In addition to the original frequency spectrum and images 25a and 25b. FIGURE 5 also illustrates an exemplary 
frequency response 29 for a suitably designed digital FIR filter having the general frequency response function, H(y), 

20 discussed above. The frequency response 29 has a shape that is defined by: a main lobe 40; and a series of side 
lobes 36a and 36b. As is well known in the filtering art, the main lobe 40 comprises a pass band 30 and a pair of 
transition bands 32a and 32b located on either side of the pass band 30. Transitions between the pass band 30 and 
the transition bands 32a and 32b are defined by points 62a and 62b. respectively The series of side lobes 36a and 
36b form a pair of stop bands 34a and 34b, located adjacent to transition bands 32a and 32b, respectively Transitions 

25 between the stop band 34a and the transition band 32a, and between the stop band 34b and the transition band 32b 
are defined by points 60a and 60b, respectively The pass band 30, transition bands 32a and 32b, and stop bands 34a 
and 34b are indicated by brackets in FIGURE 5. For ease of illustration, the frequency response 29 is centered at the 
zero frequency point. 



As is well known in the filtering art, the shape of the frequency response 29 is determined, In part, by the over- 



30 sampling factor, m. More specifically, the relationship between the transition points: 60a and 62a; and, 60b and 62b Is 
fixed and Is equal to: 2m - 1 ; where m > 1 . For example, an oversampling factor of two (I.e., m = 2) causes points 60a 
and 60b to occur at three times the frequency of points 62a and 62b. That is, the stop bands 34a and 34b begin at 
three times the cutoff frequencies of the pass band 30. In the above example (i.e. , m = 2), the pass band 30 lies between 
±fs/4 and the stop bands begin at ±3ts/4S Since the bandwidth of the original signal frequency spectrum 24 is equal 

35 to fg/m, in the above example (i.e. , m = 2), the bandwidth is equal to iJ2 and lies between ± fs/4. As a result, the original 
signal frequency spectrum 24 lies within the pass band 30 of the digital FIR filter and the images 25a and 25b lie within 
the stop bands 34a and 34b. 

If the digital FIR filter were a perfect filter, the stop bands 34a and 34b would be perfectly attenuating, and the 
frequency information passed by the digital FIR filter would simply be the frequency spectrum 24. However, because 

40 real filters are not perfect, the stop bands 34a and 34b have imperfect attenuation capabilities. The series of side lobes 
36a and 36b depict the imperfect attenuation of the stop bands 34a and 34b. The peak imperfect attenuation level is 
illustrated in FIGURE 5 by the horizontal dashed lines 44a and 44b. Because the stop bands 34a and 34b are not 
perfectly attenuating, portions of the spectral images 25a and 25b that lie within the stop bands 34a and 34b and are 
below the attenuation levels 44a and 44b are passed by the digital FIR filter. These portions of the images 25a and 

45 25b passed by the digital FIR filter produce aliasing errors in the frequency Information passed by the filter. Thus, the 
passed frequency information cannot be used to accurately replicate the original signal frequency spectrum 24. As will 
be better understood from the following discussion, the general shape of the frequency response 29 can be altered 
through appropriate design of the digital FIR filter. For example, the amplitudes of the series of side lobes 36a and 36b 
can be reduced, thereby increasing the attenuation levels 44a and 44b of the stop bands 34a and 34b. By increasing 

50 the stop band attenuation, a smaller portion of the spectral images 25a and 25b will be passed by the filter and the 
aliasing errors associated with the original signal frequency spectrum 24 passed by the digital FIR filter will be reduced. 
As a result, original signal frequency spectrum informatbn contained in the interpolated values, V|y. will contain less 
image information, so that the frequency domain information contained In the V|v values can be used to more accurately 
replicate the original signal frequency spectrum. 

55 Another source of errors that can distort the original signal frequency spectrum 24 passed by the digital FIR filter 

is the shape of the pass band 30. In an ideal filter, the pass band 30 would be perfectly flat with respect to frequency 
However, in practice, filters are not ideal and the pass band 30 is not perfectly flat. The digital FIR filter illustrated in 
FIGURE 5 has a pass band ripple 38 that may introduce amplitude errors into the frequency spectrum 24 that Is passed 
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by the filter. As will be better understood from the following discussion, the pass band flatness can be improved through 
appropriate filter design by reducing the pass band ripple 38. Reducing the pass band ripple 38 reduces the amplitude 
errors introduced into the frequency spectrum 24 by the ripple. Thus, In accordance with the invention, it is desirable 
to design a digital FIR filter that has a suitably flat pass band 30 and suitably attenuating stop bands 34a and 34b. An 

5 exemplary method for designing such a digital FIR filter suitable for use in the present invention is discussed next. 

In the preferred embodiment of the invention, as was noted above, the digital FIR filter is actually implemented in 
software form within the signal processor 10. As is well known in the filtering art, the frequency response function, H 
(y), of such a digital FIR filter depicted by the exemplary frequency response 29 in FIGURE 5, can be represented by 
one or more mathematical equations characterized by poles and zeros in the complex-frequency plane (hereinafter 

10 referred to as the frequency plane). As will be better understood from the following discussion, relocating one or more 
pairs of zeros within the frequency plane is one method of modifying the frequency shape 29. More specifically, the 
amplitudes of the side lobes 36a and 36b may be adjusted (i.e., reduced) by moving zeros along the real axis in the 
frequency plane, and the amplitude of the pass band ripple 38 may be adjusted (i.e., reduced) by moving zeros along 
the imaginary axis of the frequency plane. 

IS As is well known in the filtering art, the frequency response function, H(y), of a filter, such as a digital FIR filter, is 

related to the continuous Impulse response function, h(x), of the filter. Specifically, the h(x) function is the inverse 
Fourier transform of the H(y) function. Thus, just as with the H(y) function, the h(x) function can be represented by one 
or more mathematical equations that are the inverse Fourier transform of the equations representing the H(y) function. 
As a result, changes made to the H(y) function which alter the frequency response 29 will cause corresponding changes 

20 in the continuous impulse response function, h(x). As will be presented in the following discussion, each pair of zeros 
that is moved about the frequency plane adds a cosine term to the equations representing the impulse response 
function, h(x). Adding cosine terms to the h(x) function changes the shape of the digital FIR filter's impulse response. 
However, as will be explained more fully below, as long as the number of zeros remains the same within a suitable 
frequency band, the width of the continuous impulse response, T, remains unchanged. 

25 FIGURE 6 illustratively depicts this relationship between the H(y) and h(x) functions of the digital FIR filter. Fl GURE 

6, line A illustrates a rectangular continuous impulse response 50. Such a rectangular impulse response 50 is repre- 
sentative of the h(x) function for a perfect FIR filter. The rectangular impulse response 50 has unity height (i.e., h = 1) 
and width, T. The corresponding H(y) function for a filter having the rectangular impulse response 50 is the Fourier 
transform of the impulse response 50 and can be defined by an equation of the general form: 

30 

H(y) = ^. (1) 



The frequency response 29 represented by Equation (1) is illustrated in FIGURE 6 on line B. The zero crossings of 

35 the frequency response 29 are the zeros of Equation (1). By normalizing the frequency scale in FIGURE 6 by T the 
zeros occur at multiples of 1/T (I.e., y = ±1 . ±2, etc.). 

The frequency response illustrated in FIGURE 6, line B, can be depicted in a more conventional form by changing 
the ordinate from a linear to a logarithmic magnitude scale. This causes the negative excursions of the frequency 
response to become positive excursions while remaining in the same relative position along the frequency axis. The 

40 resulting waveform of the frequency response is illustrated in FIGURE 6 on line C and has the same general shape 
as the frequency response 29 illustrated in FIGURE 5 and discussed above. Accordingly the frequency response 
waveform illustrated in FIGURE 6, line C has a main lobe 40 and a series of side lobes 36a and 36b. For reasons of 
clarity, the pass band 30. transition bands 32a and 32b, and stop bands 34a and 34b are not illustrated in FIGURE 6. 
however, they have the same characteristics as discussed above and illustrated in FIGURE 5. 

45 FIGURE 6, tine D. is a graphical illustration of how relocating zeros along the real axis in the frequency plane can 

reduce the amplitudes of the side lobes 36a and 36b. For example, by relocating two pairs of zeros that occur at the 
normalized frequency values of ±1 and ±2 (i.e., y = ±1. ±2), the width of the main lobe 40 is increased. Moving one 
pair of zeros along the real axis causes a new pair of side lobes 46a and 46b to be inserted where the zero pair has 
been relocated, such as at y = ±3.5. for example. The reinserted side lobes 46a and 46b reduce the amplitudes of side 

50 lobes 42a and 42b, which are adjacent to the main lobe 40. In the above example, the adjacent side lobes 42a and 
42b have the highest peak value and therefore determine the attenuation levels 44a and 44b. By moving one pair of 
zeros from y = ±1 . for example, to y = ±3.5, the stop band attenuation levels 44a and 44b are increased to levels 44a* 
and 44b' and are defined by the reduced amplitudes of the adjacent side lobes 42a and 42b. 

A second pair of zeros, originally located at y = dh2 in the above example, can be relocated to a point along the 

55 Imaginary axis such as y = ± i1 .3, for example. The imaginary axis in FIGURE 6, line D is normal to the page and is 
not shown for reasons of clarity As noted above, by appropriately positioning the second pair of zeros along the im- 
aginary axis, the pass band ripple 38 (not shown) of the frequency response can be reduced. Thus, by appropriately 
relocating pairs of zeros associated with the frequency response, the pass band ripple 38 and stop band attenuation 
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levels 44a and 44b can be innproved. 

As noted above, when zeros are relocated in the frequency plane, cosine terms are added to the continuous 
impulse response function, h(x), which change the shape of the h(x) function. In the example illustrated in FIGURE 6, 
the added cosine terms cause the rectangular impulse response 50 to change to a nonrectangular impulse response 
5 52, whose shape is illustrated In FIGURE 6 on line E. The width, T. of the nonrectangular impulse response 52 remains 
unchanged because, as was noted above, the total number of zeros in the above example was not changed. 

So as to better understand the present invention, a brief discussion of the relationship between the movement of 
zero pairs in the frequency plane and the addition of cosine terms to the time domain function is presented next. As 
noted above. Equation (1 ) represents the frequency response of a perfect filter having the rectangular impulse response 
10 50 illustrated in FIGURE 6, line A. As is well known in the filtering art, a perfect filter, I.e., a rectangular response 50 
is not achievable. A more realistic filter has a nonrectangular impulse response (such as the nonrectangular Impulse 
response 52 illustrated in FIGURE 6, line E, for example). Such a nonrectangular impulse response 52 can be repre- 
sented by introducing a series of cosine terms into Equation (1 ), The corresponding frequency function. H(y), for such 
a nonrectangular impulse response 52 can be defined by an equation of the general form: 



IS 



20 



30 



3S 



40 



where: 



n is the number of cosine terms In the time domain; 
€ is an Integer; and, 
25 Is coefficient of the cosine term. 

Equation (2) represents a series of waveforms, each having the general frequency shape defined by Equation (1), 
centered at different frequencies (i.e., y - i) and having different amplitudes (i.e., a^). This is further Illustrated by writing 
Equation (2) in a slightly different form as illustrated by the following equation: 



n (-1)^ a, 

H(y) = sin ,y , (y , I) • »> 

The corresponding time domain or continuous impulse response function, h(x), is defined as the Inverse Fourier trans- 
form of Equation (2), and can be represented by an equation having the general form: 

n _ 

h(x) = R{x) I {'ly A cos 2iiix (4) 

^ where: 

R(x) is a rectangle having unit width and height. 
A^ is a coefficient of a cosine term defined as: 

50 A^ = 2(-1)%^for€;t0;and. 

Ae = for € = 0; and, 

X = ^ where T Is the width of the h(x) function. 

55 Thus, as illustrated by Equation (4), a nonrectangular impulse response 52 can be defined by a constant and a series 
of cosine terms. 

As stated above, such a nonrectangular impulse response 52 has the corresponding frequency shape defined by 
Equation (3). For reasons to be discussed below, Equation (3) can be reduced to a rational fraction having the general 
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form: 



2 2 2 2 9 2 

y-:„ ^ B,[(y^-a,)(y^^a^)...(y^-a^)] 

• Z 2 2 2~~2 (^) 



where: 



Bq is a constant; and, 

to (where € is an integer between 1 and n) is related to A« in Equation (4). The terms are the roots (i.e., zeros) 

of the numerator in Equation (5). 

As can be seen by comparing Equations (4) and (5), the number of cosine terms in the time domain (Equation (4)) is 
related to the number of new zeros in the frequency domain (Equation 5)). Moving a zero in the frequency plane is 

IS accomplished by adding an appropriate term to Equation (5) that cancels one zero term and replaces the cancelled 
zero term with a new zero term, thus relocating the cancelled zero to a new location without changing the number of 
zero terms. Adding terms to Equation (5) so as to relocate a zero causes a cosine term to be added to Equation (4) 
because, as noted above, the terms are related to the A^ terms (i.e., the cosine coefficients). 

Thus, by moving zeros of the H(y) function in the frequency plane so as to improve the stop band attenuation 44a 

20 and 44b and improve the flatness of the pass band 30, causes cosine terms to be added to the impulse response in 
the time domain. As will be readily understood from the foregoing discussion by a person of reasonable skill in this art, 
the frequency response of a digital FIR filter can be designed so as to optimize the stop band attenuation 44a' and 
44b' and pass band ripple 38 of the filter by moving zeros In the frequency plane. 

As will be better understood from the following discussion, and in accordance with the preferred embodiment of 

2S the present invention, the continuous Impulse response function, h(x), can be convolved with the data sample values, 
V3, to (eventually) produce the interpolated values, V|y, at the appropriate resample times, tp. In general, during a 
convolution operation, the h(x) function is centered at an appropriate sample time, tg, so that a tp time is encompassed 
by the finite impulse response width, T The h(x) function is multiplied by the Vs value corresponding to the centered 
{q time, thereby producing a weighted h(x) function. A convolved value, designated V^. is the particular value of the 

30 weighted h(x) function corresponding to the t^ time encompassed by the impulse response width. T. As will be better 
understood from the following discussion, a plurality of values for each t^ time are computed when a plurality of 
weighted h(x) functions overlap the Ir time. The plurality of values are summed at each tp, time to produce the 
interpolated value, V,y. The convolution steps briefly discussed above are illustrated in FIGURE 7 and discussed in 
greater detail next by way of an example. In the example which follows, a four-point digital FIR filter is used. As noted 

35 above, a four-point filter has a h(x) function whose width, T encompasses four data samples, (i.e., 4At < T < 5At). 

FIGURE 7. line A, illustrates an exemplary series of Vs values produced by a time-based sampling process. The 
Vq values occur at tg times that, preferably, are spaced apart by equal time Intervals, designated At. In the above 
example, each tg time Is assigned an integer value, designated N, starting, for example, with the first tg time (I.e., N = 
1) and increasing sequentially to the last tg time (i.e., N = Nth value). This Is illustrated in FIGURE 7. line A as the 

40 integers 1 to N above the series of Xq times. A tp time, preferably, is located between two adjacent tg points-te and t9 
or designated integers 8 and 9, for example. In this example, the four tg times located nearest to t^, designated integer 
values of seven through ten (i.e. . N = 7-1 0), are located at sample times, designated ty through t,o. and have Vg values, 
designated V7 through V^o, respectively. 

The convolution operation is Illustrated in FIGURE 7 on lines B-E. The h(x) function is represented by the nonrec- 

4S tangular impulse response shape 52. First, the impulse response 52 Is centered at ty and is multiplied by the V7 value 
to produce a convolution value, designated Vd. located at tp. Next, the impulse response 52 is centered at tg and 
multiplied by Vg to produce a convolution value, designated Vc2> located at tp. Similarly, the impulse response 52 is 
separately located at tg and t^o and multiplied by V9 and V^q to produce convolution values, designated and Vc4, 
respectively Next, the ^cv^ca values are summed to produce the V,v value at the t^ time. The convolution process 

so is repeated for subsequent t^ times until a V|v value has been computed for each tp time. 

The convolution operation illustrated in the above example can be expressed In terms of a general convolution 
equation having the following form: 



(6) 
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where: k is an integer; and, q is equal to the number of tg tinnes encompassed by the impulse response width, T, (i.e., 
q = 4 for the four-point digital FIR filter in the above example). The impulse response function, h(x), is expressed in 
terms of the tp time offset by an integer multiple of At (i.e.. Nj^At) and therefore forms an offset impulse response function 



10 Since, as was discussed above, each Xq time is assigned a sequential integer value (i.e.. 1 to N), the ts times located 
near any tp time can be defined by the following equation: 

N, = No + k (7) 

75 

where Nq is the integer multiple of At for the time closest to x^, and k is the Integer defined in Equation (6). Nq can 
be calculated from the following equation: 

No = INT (^^+ 0.5) (8) 

where Nq is rounded off to the nearest integer. 
An offset between t^ and Nq, designated t', is computed from the following equation: 

25 

t' = tR-NoAt (9) 

Solving Equation (7) for Nq and substituting this value into Equation (8) and rearranging the terms provides the following 
relationship: 

tj, - Nj^At = f - kAt. (10) 

Substituting the right-hand side of Equation (10) into the impulse response term In Equation (6) yields the following 
relationship: 



40 



(11) 



Equation (11) illustrates how the offset Impulse response function 

45 



(i.e.,h^^B_!!!^ 



(i.e.,h(-Ji-V^j ) 

SO 

of the digital FIR filter can be expressed in terms of the offset time, t' 



55 



Substituting the right side of Equation (11) into Equatbn (6) produces the following equation: 
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20 



2S 



30 



3S 



4S 



As illustrated In Equation (12), V,v can be computed from a general convolution equation by multiplying the Vg values 
(i.e., Vs(N)(At)) by the offset impulse response function 



The products are then summed to produce the V|v value at the resampling time, tpj. 
IS Generally, the value of t' is different for each tp time and. as a result, the value of 



must be computed for each tp time. There are several well-known techniques for computing 



a few of which are noted next. One technique is to store several precalculated values of the impulse response function 
(I.e., h(x)) and use these values to compute 

directly Another technique is to store only a few precalculated values of h(x) and interpolate to find 



(-11^). 



40 Additionally if a large look-up table of cosine terms is available, Equation (4) could be solved directly and the result 
could be inserted directly into Equation (6). Obviously several conventional techniques exist for computing the value of 



\ T 



and, as such, these techniques do not form a part of the present invention. 

As can be seen from examining Equation (12), the present invention provides a representation of a continuous 
approximation of the Vq signal. That is, for any t^ time desired, Equation (1 2) can be solved to provide the corresponding 

50 V|v value. As noted above, in addition to time domain information, the V|v values computed from Equation (12) contain 
frequency domain information that can be used to accurately replicate the Vq signal frequency spectrum. 

As set forth in the foregoing discussion and in accordance with the preferred embodiment of the present invention, 
a digital FIR filter is used in a novel manner to interpolate between data samples of an analog signal and produce 
interpolated values that contain accurate Vq signal frequency spectrum information. As noted above, such a digital 

55 FIR filter Is preferably implemented in software form in a program that controls a signal processor 10 suitable for 
carrying out the invention. The broad functional steps of such a program are illustrated in FIGURE 2 and discussed 
above. A more detailed flow chart of the program depicted in FIGURE 2 is illustrated in FIGURES 4a. 4b, and 4c and 
discussed next. The steps set forth in FIGURES 4a-4c are examples of a four-point digital FIR filter application, such 
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as for the four-point filter discussed in the above example. Obviously, the steps depicted in FIGURES 4a-4c can be 

readily adapted to other digital FIR filters (such as a ten-point filter, for example). 

Initially, as illustrated in FIGURE 4a, the program instructs the signal processor 10 to read and store: the V3 values; 

the Xq times; the times; the At interval; and, the data block size, Bis. The program further instructs the signal processor 
5 10 to Open a temporary file, V|y(l). for storing the computed V|y values. 

Next, the program causes the signal processor 1 0 to assign sequential integer values, N. to the tg times that have 

been read and stored. Preferably, the signal processor 10 begins with the first tg time and proceeds sequentially to the 

last tg time. That Is, for example, if there are ten ts times, they are consecutively numbered 1-10 from the first to the 

last ts time. Next, the signal processor 10 Is Instructed to read the first tp time. 
10 The next functional step performed by the signal processor 10 is to determine the set of Xq times that are closest 

to the first t,, time. As noted above, the size of the set is determined by the width of the impulse response. T. First, a 

resampling index counter, designated I, and a summation index counter, designated K, are initialized (i.e., I = 1 and K 

= 1). Next, the program instructs the signal processor 10 to determine the closest \q time to the first t^ time, i.e., Nq. 

Nq is computed from the following equation: 



20 



INT 



At 



+ 0.5 4-1 



(13) 



Next, the offset, t'. between the first resample time, t^, and the nearest tg time, (i.e., Nq • At) is computed from 
Equation (9) above. Once Nq and f have been computed, the program instructs the signal processor 10 to find the 
2S other tg times. In the set of ts times closest to t^. In the present example, I.e.. a four-point digital FIR filter application, 
the signal processor 10 will locate the four tg times closest to tp. Specifically and as illustrated in FIGURE 4b, if t' Is 
equal to or greater than half the interval between adjacent Ig times (i.e., t' > At/2), then the four closest tg times are 
determined by the following equation: 

30 

Nk = No + K-C (14) 



If t' Is less than half the At interval (i.e.. t' < At/2), then the four closest tg values are determined by the following equation: 

Nk = No+C-K (15) 

In Equations (14) and (15). K is the summation Index counter value noted above and C is an integer whose value is 
one-half the number of tg times encompassed by T (i.e., where C = 2 for the four-point filter in the above example). 

After the first N^^ point has been determined from Equation (14) or (15), the program determines whether the N,^ 
point lies within the data block (i.e., l<N,^<By and, If not, adjusts the N^ value so that it is within the data block. More 
specifically, if the integer value of the tg time lies outside the data block (i.e., N,^ > Bis), and is greater than or equal to 
one (i.e., N^^l ), the program places the tg time within the data block by subtracting the data block size from the Integer 
value of the tg time (i.e., N^ = N^ - Bis). The signal processor 10 stores this value of Nk- If the integer value of the tg 
time is outside the data block and is less than one (i.e., Nj^ < 1). then the program adds the data block size to the 
Integer value of the tg time (i.e., N^ = N^ + Bis) and stores this value of N^. If the integer value of the tg time is within 
the data block (i.e., 1<NK<By. the value of N^ Is stored. 

Next, the program determines If all of the tg times in the set have been found. Specifically, in the four-point filter 
example discussed above, if the summation index counter value is less than four (i.e., K < 4) then the counter is 
incremented (i.e., K = K + 1) and the above process is repeated. A summation index counter value greater than or 
equal to four (i.e., K ^ 4) indicates that the four tg times in the set have been found and the program proceeds to the 
next functional step, which is described next and illustrated in FIGURE 4c. 

Once the appropriate number of tg times near the t^ time have been found, the signal processor 10 convolves the 
continuous impulse response function (h(x)), of a suitably designed digital FIR filter with the Vg values of the tg times 
found in the set discussed above. As noted above, a suitable digital FIR filter is preferably implemented in software 
form and can be characterized by one or more mathematical formulae, such as Equation (2) (i.e., the frequency re- 
sponse function. H(y)) and Equation (4) (i.e., the impulse response function, h(x)). for example. The program instnicts 



11 



EP 0 355 293 B1 

the signal processor 10 to determine a particular value of the offset impulse response function, 




T )' 



defined in Equation (1 1 ) above, where k is replaced by K in the progranr^ illustrated in FIGURE 4C. As discussed above, 
the value of 



can be determined by conventional techniques. Next, the signal processor 10 performs the convolution operation by 
using the value of 



in the general convolution equation set forth in Equation (12) where, as noted above, k = K and q = 4. The resulting 
V,v value is stored as a V,v (I) value. This last step in the program performs the two functional steps illustrated in 
FIGURE 2 and discussed above. That is, the functional steps of producing the Vq values and summing them to produce 
a V,v value are performed by the signal processor 10 in the convolution operation step illustrated in FIGURE 4C. 

Once the V|v value has been computed and stored by the signal processor 10, the program increments the resa- 
mpling counter. I (i.e., I = I + 1). If the resampling counter value is less than or equal to the total number of resample 
times read by the signal processor 10 (i.e., I < Bis), then the next tp time is read and the above steps are repeated. A 
resampling counter value that exceeds the total number of resample times (i.e.. I > Bis) indicates that the V|v values 
have been computed for all t^ times read by the signal processor 10 and the program is exited. 

As can be readily appreciated from the foregoing description, the invention provides a method and apparatus tor 
preserving the frequency spectrum of a sampled analog signal using a digital FIR filter in a novel manner for Interpolating 
between data samples. As a result, the interpolated values contain frequency information that can be used to accurately 
replicate the frequency spectrum of the analog signal. While a preferred embodiment of the invention has been illus- 
trated and described herein, it is to be understood that, within the scope of the appended claims, various changes can 
be made. For example, digital FIR filters other than the four-point filter discussed in the above example can also be 
used. In fact, even non-FIR filters (i.e., filters having continuous impulse responses of infinite width) can be utilized, 
liowever, the convolution operation becomes much more complex. The frequency response and impulse response of 
the digital FIR filter can also be different than as illustrated and discussed above. Likewise, an oversampling factor 
other than two (I.e., 1 < m = 2) could also be used. While the time inten/al, At, between adjacent data samples is 
preferably the same, the invention is capable of performing its intended function even when the data samples are 
produced by a sampling process that samples the analog signal at a variable rate. Furthermore, the program merely 
illustrates one possible method of carrying out the invention and it is to be understood that other equally suitable 
programs could be utilized. 



1. A method of producing interpolated values of an analog signal which has been sampled at known sample times 
(ts). the method comprising the steps of: 

(a) storing data samples (Vs) of the anabg signal, said data samples containing frequency domain infornnation 
about the analog signal's original frequency spectrum and frequency domain information about images of the 
analog signal's original frequency spectrum; 

(b) receiving an arbitrary selection of desired resample times (tr) at which said data samples are to be inter- 
polated; and for each resample time (tr) : 

(c) choosing a set of data samples (Vs) occurlng at sample times (ts) that are located nearest the resample 
time (tr); 

(d) weighting said set of data samples (Vs) by convoluting each of said data samples (Vs) with a continuous 





Claims 
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nonrectangular impulse response function h(x) ol a filter having a predetermined frequency response H(y) 
comprising a substantially flat passband and suitably attenuating stopbands, said weighting producing a plu- 
rality of convolution values at each of said arbitrary selection of resampte times (tr); and, 
(e) summing said plurality of convolution values at each of said resample times and producing an interpolated 
value at each of said resample times (tr). 

Method according to claim 1, wherein the size of the set of sample values (Vs) is determined by the width of the 
impulse response of the filter. 

The method according to one of the preceding claims, wherein the analog signal is sampled at a sampling fre- 
quency, fs, with an oversampiing factor, m, and the predetemnined frequency response includes a substantially 
flat passband for a band ot frequencies between - fs/2m and fs/2m, and suitably attenuating stopbands for fre- 
quencies above ts-^^J^ and below -fs-^^^j^. 

Method according to one of claims 1 to 3, wherein a frequency function for the continuous nonrectangular impulse 
response has a general form as follows: 




An apparatus for producing interpolated values of an analog signal which has been sampled at known sampling 
times (ts), the apparatus comprising: 

(a) a receiver (18) for receiving data samples (Vs) of the analog signal, said data samples contain frequency 
domain information about the analog signal's original frequency spectrum and frequency domain information 
about images of the analog signal's original frequency spectrum, and for receiving an arbitrary selection of 
resample times (tr) at which said data samples (Vs) are to be interpolated; and, 

(b) interpolating means (20, 21 , 22, 23) coupled to said receiver (18) for receiving said data samples (Vs) and 
said selection of resample times (tr), said interpolating means Interpolating between said data samples to 
produce interpolated values (Viv) at said resample times (tr), wherein said interpolating means includes 

(b1) a filter (21) producing a continuous nonrectangular impulse response h(x) and an associated fre- 
quency response (Hy) comprising a substantially flat pass-band and suitably attenuating stopbands; 
(b2) a data locator (20) coupled to said receiver ( 1 6) for receiving said data samples (Vs) and said resample 
times (tr) and locating a set of said data samples (Vs) for each resample time (tr) that is located nearest 
each of said resample times (tr); 

(b3) convolving means (22), coupled to said data filter (21) and said data locator (20), for weighting the 
set of the data samples (Vs) with the continuous nonrectangular impulse response of the filter (21) for 
each of said resample times (tr) to produce a plurality of convolution values at the respective resample 
times (tr); and 

(b4) a summer (23) coupled to the convolving means (22) which sums the plurality of convolution values 
at each of the resample times to produce one of the interpolated values (Viv) at each of the resample 
times (tr). 

The apparatus according to claim 5 wherein a frequency function for the continuous nonrectangular impulse re- 
sponse has a general form as follows: 

The apparatus according to claim 5 or 6, wherein said frequency response of said filter (21 ) comprises: 

(a) a pass band for passing frequency domain information about the frequency spectrum of the analog signal 
and reducing contribution from pass band errors contained in said interpolated values; and, 

(b) stop bands for reducing contribution from Images of the frequency spectrum of the analog signal contained 
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in said interpolated values. 
8. The apparatus according to one of claims 5 to 7, wherein said filter (21 ) is a digital filter. 
s 9. The apparatus according to one of claims 5 to 8, wherein said filter (21 ) is a digital finite impulse response filter. 
10. The apparatus according to one of claims 5 to 9, wherein said filter (21) is implemented in software form. 

10 Patentanspruche 

1. Verfahren zum Erzeugen interpoiierter Werte eines analogen Signals, das bel bekannten Abtastzeiten (ts) abge- 
tastet wurde, mit folgenden Verfahrensschritten: 

'5 (a) Spelchern der Datenabtastwerte (Vs) des analogen Signals, wobei die Datenabtastwerte Frequenzbe- 

reichsinformation uberdas ursprungiiche Frequenzspektrum des analogen Signals und Frequenzbereichsln- 
formation uber Bilder des ursprunglichen Frequenzspektrums des analogen Signals enthalten; 

(b) Empfangen einer willkOrllchen Auswahl envunschter erneuter Abtastzeiten (tr), zu denen die Datenabtast- 
werte interpoliert werden sollen; und 

20 fur jede erneute Abtastzeit (tr): 

(c) Wahlen einer Gruppe aus Datenabtastwerten (Vs), welche bei den Abtastzeiten (ts) auftreten, die am 
nachsten bei der erneuten Abtastzeit liegen; 

(d) Gewichten der Gruppe der Datenabtastwerte (Vs) durch Fallen jedes Datenabtastwertes (Vs) mit einer 
kontinuierlichen nicht rechteckigen Impulsantwortfunktion h(x) eines Filters mit einem vorgegebenen Fre- 

25 quenzgang H(y), welcher einen Im wesentlichen flachen Durchla3bereich und geeignet dampfende Sperrbe- 

reiche aufweist, wobei das Gewichten mehrerer Faltungswerte bei jeder der willkurlich ausgewahlten erneuten 

Abtastzeiten (tr) erzeugt; und 

(e) Summieren der mehreren Faltungswerte bei jeder der erneuten Abtastzeiten und Erzeugen eines interpo- 
lierten Wertes bei jeder der erneuten Abtastzeiten (tr). 

30 

2. Verfahren nach Anspruch 1 , bei dem die Gr63e der Gruppen der Abtastwerte (Vs) durch die Breite der Impulsant- 
wort des Filters bestimmt wird. 

3. Verfahren nach einem der vorangehenden Anspruche. bei dem das analoge Signal mit einer Abtastfrequenz, fs, 
3S mit einem Oberabtastungsfaktor, m, abgetastet wird, und die vorgegebene Frequenzantwort einen im wesentlichen 

flachen Durchlaf3bereich fur ein Frequenzband zwischen -fs/2m und fs/2m sowie geeignet dampfende Sperrbe- 
relche fur Frequenzen uber fs^^i^ und unter -ts^^J^ aufweist. 

4. Verfahren nach einem der AnsprOche 1 bis 3, bei dem eine Frequenzfunktion fOr die kontinuierliche nicht recht- 
40 eckige Impulsantwort die folgende allgemeine Form hat: 



45 



5. Vorrlchtung zum Erzeugen interpoiierter Werte eines analogen Signals, welches bel bekannten Abtastzeiten (ts) 
abgetastet wurde, mit folgenden Merkmalen: 

(a) einem Empfanger (18) zum Empfangen von Datenabtastwerten (Vs) des analogen Signals, wobei die 
Datenabtastwerte Frequenzbereichsintonmation uber das ursprungiiche Frequenzspektrum des analogen Si- 
gnals und Frequenzbereichsintormation uber Bilder des ursprunglichen Frequenzspektrums des analogen 
Signals enthalten, und zum Empfangen einer willkOrllchen Auswahl erneuter Abtastzeiten (tr), bei welchen die 
Datenabtastwerte (vs) interpoliert werden sollen; 
55 (b) Interpolationsmittel (20, 21 , 22, 23), welche mit dem Empfanger (18) verbunden sind, zum Empfangen der 

Datenabtastwerte (Vs) und der erneuten Abtastzeiten (tr) wobei die Interpolationsmittel zwischen den Daten- 
abtastwerten interpolieren, um interpolierte Werte (Viv) bei den erneuten Abtastzeiten (tr) zu erzeugen. wobei 
die Interpolationsmittel folgende Merkmaie aufweisen: 
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(b1) ein Filter (21), welches eine kontinuierliche nicht rechteckige Impulsantwort h(x) und einen zugeho- 
rigen Frequenzgang H(y) mit einem Im wesentlichen flachen DurchlaBbereich und geeignet dampfenden 
Sperrbereichen erzeugt; 

(b2) einen Datenlokalisierer (20), der nnit dam Empfanger (18) verbunden ist, zum Emplang der Daten- 
abtastwerte (Vs) und der erneuten Abtastzeiten (tr) und zum Lokalisieren einer Gruppe aus Datenabtast- 
werten (Vs) fur jede erneute Abtastzeit (tr), welche am nachsten bei jeder der erneuten Abtastzeiten (tr) 
liegt, 

(b3) Faltungsmittel (22). welche mit dem Datenfilter (21) und dem Datenlokalisierer (20) verbunden sind, 
zum Gewichten der Gruppe der Datenabtastwerte (Vs) mit der kontinuierlichen nicht rechteckigen Impuls- 
antwort des Filters (21 ) fur jede der erneuten Abtastzeiten (tr). um mehrere Faltungswerte bei den jewei- 
ligen erneuten Abtastzeiten (tr) zu erzeugen, 

(b4) einen Summierer (28), der mit den Faltungsmittein (22) verbunden ist und die mehreren Faltungswerte 
bei jeder der erneuten Abtastzeiten summiert, um bei jeder der erneuten Abtastzeiten (tr) einen der inter- 
polierten Werte (Viv) zu erzeugen. 

6. Vorrichtung nach Anspruch 5, bei der die Frequenzfunktion fur die kontinuierliche nicht rechteckige Impulsantwort 
die folgende allgemeine Form hat: 



7. Vorrichtung nach Anspruch 5 Oder 6, bei der der Frequenzgang des Filters (21 ) folgende Merkmale aufweist: 

(a) einen Durchla3bereich zum Durchlassen von Frequenzbereichsinformationen Oberdas Frequenzspektrum 
des analogen Signals und zum Reduzieren von Beitragen aus DurchiaBbereichsfehlern, welche in den inter- 
polierten Werten enthalten sind; und 

(b) Sperrbereiche zum Reduzieren von Beitragen aus Bildern des Frequenzspektrums des analogen Signals, 
welche In den interpolierten Werten enthalten sind. 

8. Vorrichtung nach einem der Anspruche 5 bis 7. bei der das Filter (21 ) ein digitales Filter ist. 

9. Vorrichtung nach einem der Anspruche 5 bis 8, bei der das Filter (21 ) ein digitales FIR-Filter ist. 

10. Vorrichtung nach einem der Anspruche 5 bis 9, bei der das Filter (21) in Softwareform realisiert ist. 



Revendicatlons 

1 . Proced6 pour produire des valeurs lnterpol6es d'un signal analogique qui a 6t6 echantillonn6 k des temps d'6chan- 
tlllonnage conn us (ts), le proc6d6 comprenant les Stapes consistant k: 

(a) mdmoriser des dchantiltons de donnSes (Vs) du signal analogique, lesdits Schantillons de donnSes con- 
tenant une information de domatne fr6quentiel concernant le spectre de frequence d'orlgine du signal analo- 
gique et une information de domaine f rSquentiel concernant des images du spectre de frequence d'origine du 
signal analogique; 

(b) recevoir une selection arbitraire de temps de reechantillonnage souhaites (tr) auxquels lesdits echantillons 
de donn6es doivent §tre interpoISs; et 

pour chaque temps de reechantillonnage (tr): 

(c) choisir un ensemble d'dchantillons de donnees (Vs) apparaissant k des temps d'dchantillonnage (ts) qui 
sont situds les plus proches du temps de reechantillonnage (tr); 

(d) ponderer ledit ensemble d'echantillons de donnees (Vs) en convoluant chacun desdits echantillons de 
donnees (Vs) avec une fonction de reponse impulsionnelle non rectangulaire continue h(x) d'un filtre ayant 
une reponse en frequence predeterminee H(y) comprenant une bande passante sensiblement plate et des 
bandes d'arret avec une attenuation appropri6e, ladite ponderation produisant une plurality de valeurs de 
convolution k chacun des temps de ladite selection arbitraire de temps de reechantillonnage (tr); et 

(e) additionner ladite plurality de valeurs de convolution k chacun desdits temps de reechantillonnage et pro- 
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duire une valeur interpol6e ^ chacun desdits temps de r66chantillonnage (tr). 

Precede selon la revendication 1 , dans lequel la taille de ['ensemble de valeurs d'echantillons (Vs) est ddtemiln^e 
par la largeur de la reponse impulsionnelle du flltre. 

ProcSdd selon I'unequelconque des revendications pr^c^dentes, dans lequel le signal analogique est dchantiilon- 
n6 ^ une frequence d'echantillonnage fs avec un lacteur de surechantillonnaga m, et la r6ponse en frequence 
pr6d6termin6e inclut une bande passante sensiblement plate pour une bande de frequence entre -y2m et y2m. 
et des bandes d'arrdt avec une attenuation appropriSe pour les frequences sup^rleures ^ fs(2m-1 )/2m et inf Srieures 
^-fs(2b-1)/2m. 

Precede sebn Tune quelconques des revendications 1^3, dans lequel une fonction de frequence pour la reponse 
Impulsionnelle non rectangulaire continue a la forme gdn^rale suivante: 



n 

H(y) = 2 aJsin7r(y-L)]/[?r(y-L)]. 
L=-n 



Dispositif pour produire des valeurs interpoiees d'un signal analogique qui a dchantillonne h des temps d'echan- 
tillonnage connus (ts), le dispositif comprenant: 

(a) un recepteur (1 8) pour recevoir des echantillons de donnees (Vs) du signal analogique, lesdits echantillons 
de donnees contenant une Information de domaine frequentiel concernant le spectre de frequence d'origine 
du signal analogique et une information de domaine frequentiel concernant des images du spectre de fre- 
quence d'origine du signal analogique, et pour recevoir une selection arbitraire de temps de reechantillonnage 
(tr) auxquels lesdits echantillons de donnees (Vs) doivent etre interpoles; et 

(b) des moyens d'interpolation (20. 21, 22, 23) coupl6s audit recepteur (18) pour recevoir lesdits echantillons 
de donnees (Vs) et ladite selection de temps de reechantillonnage (tr), lesdits moyens d'interpolation realisant 
une interpolation entre lesdits echantillons de donnees pour produire des valeurs interpoiees (V|v) auxdits 
temps de reechantillonnage (tr), dans lequel lesdits moyens d'interpolation incluent 

(b1) un filtre (21 ) produisant une reponse impulsionnelle non rectangulaire continue h(x) et une reponse 
en frequence (Hy) associee comprenant une bande passante sensiblement plate et des bandes d'arrfit 
avec une attenuation appropriee; 

(b2) un circuit de localisation de donnees (20) couple audit recepteur (18) pour recevoir lesdits echantillons 
de donnees (Vs) et lesdits temps de reechantillonnage (tr) et localiser un ensemble desdits echantillons 
de donnees (Vs) pour chaque temps de reechantillonnage (tr) qui est srtue le plus proche de chacun 
desdits temps de reechantillonnage (tr); 

(bS) des moyens de convolution (22), couples audit filtre de donnees (21) et audit circuit de localisation 
de donnees (20), pour ponderer I'ensemble des echantillons de donnees (Vs) avec la reponse impulsion- 
nelle non rectangulaire continue du filtre (21 ) pour chacun desdits temps de reechantillonnage (tr) afin de 
produire une plurallte de valeurs de convolution aux temps de reechantillonnage (tr) respectlfs; et 
(b4) un additionneur (23) couple aux moyens de convolution (22) qui additionnent la pluralite de valeurs 
de convolution. ^ chacun des temps de reechantillonnage, pour produire une des valeurs interpoiees (V|v) 
k chacun des temps de reechantillonnage (tr). 

Dispositif selon la revendication 5, dans lequel une fonction de frequence pour la reponse impulsionnelle non 
rectangulaire continue a la forme generale suivante: 



n 

H(y) = 2 QrJsin7r(y-L)]/[7r(y-L)], 



Dispositif selon la revendication 5 ou 6. dans lequel ladite reponse en frequence dudit filtre (21) comprend: 
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(a) une bande passante pour transmettre une information de domalne Ir^quentiel concernant le spectre de 
frequence du signal analogique et reduire une contribution provenant d'erreurs de bande passante contenues 
dans lesdites valeurs interpoiees; et 

(b) des bandes d'arret pour reduire une contribution provenant d'images du spectre de frequence du signal 
analogique contenues dans lesdites valeurs interpoiees. 

8. Dispositif selon I'une quelconque des revendlcations 5^7, dans lequel ledit filtre (21) est un filtre numdrique. 

9. Dispositif selon I'une quelconque des revendlcations 5^8. dans lequel ledit filtre (21) est un filtre numdrique k 
r^ponse impulsionnelle finie. 

10. Dispositif selon I'une quelconque des revendications 5^9, dans lequel ledit filtre (21) est implants sous la forme 
d'un logiciel. 
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